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I. INTRODUCTION. 

Entangled states have different remarkable applications and 
among them are quantum cryptography [1, 2], superdense 
coding [3, 4], teleportation [5, 6] and the potential speedup 
of quantum algorithms [7, 8, 9]. The entanglement of bipar- 
tite systems is well-understood [10, 11, 12, 13], while the en- 
tanglement of multipartite systems offers a real challenge to 
physicists. In contrast to bipartite setting, there is no unique 
treatment of the maximally entangled states for multipartite 
systems. In this reason it is highly difficult to formulate a the- 
ory of multipartite entanglement. Another point which makes 
difficult to understand the entanglement for the multi-qubit 
systems is mainly due to the fact that the analytic expressions 
for the various entanglement measures is extremely hard to 
derive. 

We consider pure three qubit systems [14, 15, 16, 17], al- 
though the entanglement of mixed states attracts a consider- 
able attention. For example, in recent experiment [18] the tan- 
gle for general mixed states was evaluated, which has never 
been done before. Three-qubit system is important in the 
sense that it is the simplest system which gives a non-trivial 
effect in the entanglement. Thus, we should understand the 
general properties of the entanglement in this system as much 
as possible to go further more complicated higher qubit sys- 
tem. The three-qubit system can be entangled in two inequiv- 
alent ways GHZ [19] and W, and neither form can be trans- 
formed into the other with any probability of success [20]. 
This picture is complete: any fully entangled state is SLOCC 
equivalent to either GHZ or W. 

Only very few analytical results for tripartite entanglement 
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have been obtained so far [21] and we need more light on 
the subject. This is our main objective and we choose geo- 
metric measure of entanglement E g [22, 23, 24, 25]. It is an 
axiomatic measure [22, 26, 27, 28], is connected with other 
measures [29, 30] and has an operational treatment. Namely, 
for the case of pure states it is closely related to the Grove- 
rian measure of entanglement [31] and the latter is associated 
with the success probability of Graver's search algorithm [32] 
when a given state is used as the initial state. 

Geometric measure depends on entanglement eigenvalue 
A^ ax and is given by formula E g (ip) = 1 — A^ nax . For pure 
states the entanglement eigenvalue is equal to the maximal 
overlap of a given state with any complete product state. The 
maximization over product states gives nonlinear eigenprob- 
lem [25] which, except rare cases, does not allow the complete 
analytical solutions. 

Recently the idea was suggested that nonlinear eigenprob- 
lem can be reduced to the linear eigenproblem for the case of 
three qubit pure states [33]. The idea is based on theorem stat- 
ing that any reduced (n — 1) -qubit state uniquely determines 
the geometric measure of the original n-qubit pure state. This 
means that two qubit mixed states can be used to calculate the 
geometric measure of three qubit pure states and this will be 
fully addressed in this work. 

The method gives two algebraic equations of degree six 
defining the geometric measure of entanglement. Thus the 
difficult problem of geometric measure calculation is reduced 
to the algebraic equation root finding. Equations contain valu- 
able information, are good bases for the numerical calcula- 
tions and may test numerical calculations based on other nu- 
merical techniques [9]. 

Furthermore, the method allows to find the nearest separa- 
ble states for three qubit states of most interest and get ana- 
lytic expressions for their geometric measures. It turn out that 
highly entangled states have their own feature. Each highly 
entangled state has a vicinity with no product state and all 
nearest product states are on the boundary of the vicinity and 
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form an one-parametric set. 

In Section II we derive algebraic equations defining the ge- 
ometric entanglement measure of pure three qubit states and 
present the general solution. In Section III we examine W- 
type states and deduce analytic expression for their geometric 
measures. States symmetric under permutation of two qubits 
are considered in Section IV, where the overlap of the state 
functions with the product states are maximized directly. In 
last Section V we make concluding remarks. 

II. ALGEBRAIC EQUATIONS. 

We consider three qubits A,B,C with state function \tp). 
The entanglement eigenvalue is given by 



n = tr(p A cr), r 2 = ti{p B a), gij = tv{p AB o l <g> aj) (6) 

and <7,'s are Pauli matrices. The matrix is not necessarily 
to be symmetric but must has only real entries. The maxi- 
mization gives a pair of equations 

ri + gs 2 = Aisi, r 2 + g T s 1 = X 2 s 2 , (7) 

where Lagrange multipliers Ai and A 2 are enforcing unit na- 
ture of the Bloch vectors. The solution of Eq.(7) is 

si = (AiA 2 l -gg T y l (A 2 n + firr 2 ) , (8a) 



J^max = max \(q 1 q 2 q 3 \*p)\ 



(1) 



and the maximization runs over all normalized complete prod- 
uct states Iq 1 ) ® \q 2 ) ® \q 3 ). Superscripts label single qubit 
states and spin indices are omitted for simplicity. Since in 
the following we will use density matrices rather than state 
functions, our first aim is to rewrite Eq.(l) in terms of density 
matrices. Let us denote by p ABC = \ip)(ijj\ the density matrix 
of the three-qubit state and by g k = \q k ) (q k \ the density ma- 
trices of the single qubit states. The equation for the square of 
the entanglement eigenvalue takes the form 
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An important equality 



maxtr(p ABC g 1 



max tr (p ABC g 1 



(2) 



? 3 ) 



tr(/ s V 



g 2 ®! 3 ) (3) 



was derived in [33] where 1 is a unit matrix 
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meaning. The matrix tr(p ABC g 1 



It has a clear 
g") is 2 (g) 2 hermitian 
matrix and has two eigenvalues. One of eigenvalues is al- 
ways zero and another is always positive and therefore the 
maximization of the matrix simply takes the nonzero eigen- 
value. Note that its minimization gives zero as the minimiza- 
tion takes the zero eigenvalue. 

We use Eq.(3) to reexpress the entanglement eigenvalue by 
reduced density matrix p AB of qubits A and B in a form 



s 2 = (AiA 2 l - g T g) 1 (Air 2 + g T r x ) . 



(8b) 



Now, the only unknowns are Lagrange multipliers, which 
should be determined by equations 



|sif = l, |S2| 2 = 1. 



(9) 



In general, Eq.(9) give two algebraic equations of de- 
gree six. However, the solution (8) is valid if Eq.(7) sup- 
ports a unique solution and this is by no means always the 
case. If the solution of Eq.(7) contains a free parameter, then 
det(AiA 2 l — gg T ) = and, as a result, Eq.(8) cannot not 
applicable. The example presented in Section III will demon- 
strate this situation. 

In order to test Eq.(8) let us consider an arbitrary superpo- 
sition of W 



|W0 = -^(|100> + |010> + |001>) do) 



and flipped W 



\W)=^= (|011> + |101> + |110» (11) 



states, i.e. the state 



A 



< (^)=maxtr (p AB g 1 

g 1 g 2 



(4) 



We denote by Si and s 2 the unit Bloch vectors of the den- 
sity matrices g 1 and g 2 respectively and adopt the usual sum- 
mation convention on repeated indices i and j. Then 



\ip) = cos 6\W) + sin 6\W). 
Straightforward calculation yields 



n = r 2 = — (2 sin 29i + cos 29n) . 



(12) 



(13a) 



Amax = T , ma , x (1 + *i ' r i + s 2 ■ r 2 + g tj s u s 2j ) , (5) 

4 s? = si = l 



where 



1 



'2 



9 = 7 2 | , 
6 \0 0-1, 



(13b) 
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where unit vectors i and n are aligned with the axes x and z, 
respectively. Both vectors i and n are eigenvectors of matri- 
ces g and g T . Therefore Si and s 2 are linear combinations 
of i and n. Also from ri = r 2 and g = g T it follows that 
Si = s 2 and Ai = A 2 . Then Eq.(8) for general solution give 



Si = cosan + sin am, s 2 — cos f3n + sin /3 m. (21) 
Then Eq.(7) reduces to the following four equations 



where 



Si = s 2 = sin 2ip i + cos 2ip n 



. „ 2sin26» „ cos26> 

sm2ip = — — , cos2(,5 = 



3A-2 



3A + 1 



(14) 



(15) 



The elimination of the Lagrange multiplier A from Eq.(15) 
gives 



r\ — r 3 cos (3 — Ai cos a, r 2 — r 3 cos a = A 2 cos (3, (22a) 



uj sin (3 = Ai sin a, w sin a = A 2 sin (3, 



(22b) 



which are used to solve the four unknown constants Ai, A 2 , a 
and (3. Eq.(22b) impose either 



AiA 2 - uj 2 = 



or 



(23) 



3 sin 2ip cos 2ip — cos 26 sin 2ip — 2 sin 20 cos 2cp. (16) 

Let us denote by t = tan (p. After the separation of the 
irrelevant root t = — tan 9, Eq.(16) takes the form 



sintft 3 + 2cos9t 2 - 2sin9t - cosfl = 0. (17) 

This equation exactly coincides with that derived in [25]. 
Since a detailed analysis was given in Ref.[25], we do not 
want to repeat the same calculation here. Instead we would 
like to consider the three-qubit states that allow the analytic 
expressions for the geometric entanglement measure by mak- 
ing use of Eq.(7). 



III. W-TYPE STATES. 



Consider W-type state 



sin a sin f3 = 0. 



(24) 



First consider the case n > 0, r 2 > 0, r 3 > and coef- 
ficients a, b, c form an acute triangle. Eq.(24) does not give 
a true maximum and this can be understood as follows. If 
both vectors Si and s 2 are aligned with the axis z, then the 
last term in Eq.(5) is negative. If vectors Si and s 2 are an- 
tiparallel, then one of scalar products in Eq.(5) is negative. In 
this reason A^ ax cannot be maximal. Then Eq.(23) gives true 
maximum and we have to choose positive values for Ai and 
A 2 to get maximum. 

First we use Eq.(22a) to connect the angles a and (3 with 
the Lagrange multipliers Ai and A 2 



cos a = 



A 2 ri - r 2 r 3 



cos (3 = 



A x r 2 - rir 3 



(25) 



Then Eq.(22b) and (23) give the following expressions for 
Lagrange multipliers Ai and A 2 



= a|100) + 6|010) + c|001), a 2 + b 2 + c 2 = 1. (18) 

Without loss of generality we consider only the case of posi- 
tive parameters a,b,c. Direct calculation yields 



n = r t n, r 2 = r 2 n, 




(19) 



where 

ri = b 2 + c 2 -a 2 , r 2 = a 2 + c 2 -b 2 , r 3 = a 2 + b 2 -c 2 (20) 

and lo = 2ab. The unit vector n is aligned with the axis z. 
Any vector perpendicular to n is an eigenvector of g with 
eigenvalue u. Then from Eq.(7) it follows that the compo- 
nents of vectors S\ and s 2 perpendicular to n are collinear. 
We denote by m the unit vector along that direction and pa- 
rameterize vectors Si and s 2 as follows 



Al = LU 



A 2 = uj 



OJ 



1/2 



uj z +ri 



2\ 1/2 



uj 2 + r\ - r\ 



(26a) 



(26b) 



Eq.(7) allows to write a shorter expression for the entangle- 
ment eigenvalue 



A max = ^(l + A 2 +ri cosa). 



(27) 



Now we insert the values of A 2 and cos a into Eq.(27) and 
obtain 



4A 2 



<\/{uj 2 + r\- rl)(uj 2 + r 2 - rf) - r x r 2 r 3 



uj^ — rS 



(28) 
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The denominator in above expression is multiple of the area 
S of the triangle a, b, c 

lo 2 - r\ = 16S 2 . (29) 
A little algebra yields for the numerator 

uj\j {uj 2 +r\- r§) + (uj 2 + r\ - r\) - nr 2 r 3 (30) 
= I6a 2 b 2 c 2 -uj 2 + r\. 

Combining together the numerator and denominator, we 
obtain the final expression for the entanglement eigenvalue 

AL X = 4R 2 , (31) 

where R is the circumradius of the triangle a,b,c. Entangle- 
ment value is minimal when triangle is regular, i.e. for W-state 
andA 2 nax (W)=4/9[25,34]. 

Now consider the case r 3 < 0. Since r 3 + n = 2b 2 > 0, 
we have n > and similarly r2 > 0. Eq.(24) gives true 
maximum in this case and both vectors are aligned with the 
axis z 



si = s 2 = n (32) 
resulting in A 2 ^ = c 2 . In view of symmetry 

A max = max(a 2 , b 2 , c 2 ), max(a 2 , b 2 , c 2 ) > i. (33) 

Since the matrix g and vectors n and r 2 are invariant under 
rotations around axis z the same properties must have Bloch 
vectors s\ and s 2 . There are two possibilities: 

i)Bloch vectors are unique and aligned with the axis z. The 
solution given by Eq.(32) corresponds to this situation and 
the resulting entanglement eigenvalue Eq.(33) satisfies the in- 
equality 



have same entanglement value and it seems to imply some- 
thing interesting. GHZ state can be used for teleportation and 
superdense coding, but W-state cannot be. However, the W- 
type state with right triangle coefficients can be used for tele- 
portation and superdense coding [35]. In other words, both 
type of states can be applied provided they have the required 
entanglement eigenvalue A 2 ^ = 1/2. 

IV. SYMMETRIC STATES. 

Now let us consider the state which is symmetric under per- 
mutation of qubits A and B and contains three real indepen- 
dent parameters 

\ip) = a|000> + + c|001) + d|110), (36) 

where a 2 + b 2 + c 2 + d 2 = 1. According to Generalized 
Schmidt Decomposition [14] the states with different sets of 
parameters are local-unitary(LU) inequivalent. The relevant 
quantities are 

(uj 0\ 

n = r 2 = rn, g = -w , (37) 

V° V 

where 

r = a 2 + c 2 - b 2 - d 2 , uj = 2ad+ 2bc (38) 

and the unit vector n again is aligned with the axis z. 
All three terms in the l.h.s. of Eq.(5) are bounded above: 

• si ■ ri < |r|, 

• s 2 -r 2 < \r\, 

• and owing to inequality \uj\ < 1, gij sus 2 j < 1. 

Quite surprisingly all upper limits are reached simultane- 
ously at 



\ < AL X < 1- (34) 

ii)Bloch vectors have nonzero components in xy plane and 
the solution is not unique. Eq.(21) corresponds to this situa- 
tion and contains a free parameter. The free parameter is the 
angle defining the direction of the vector m in the xy plane. 
Then Eq.(31) gives the entanglement eigenvalue in highly en- 
tangled region 

\<^<\- (35) 

Eq.(31) and (33) have joint curves when parameters a, b, c 
form a right triangle and give A 2 lax = 1/2. The GHZ states 



s\= s 2 = Sign(r)n, (39) 

which results in 



A max = 2 (! + M) ■ ( 40 ) 

This expression has a clear meaning. To understand it we 
parameterize the state as 

|V> = «!i|00gi) + fc2|ll< & ), (41) 

where qi and q 2 are arbitrary single normalized qubit states 
and positive parameters k\ and k 2 satisfy k\ + k\ = 1. Then 
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AL x = max(fc2 ; fc2 )j (42) 

i.e. the maximization takes a larger coefficient in Eq.(41). In 
bipartite case the maximization takes the largest coefficient 
in Schmidt decomposition [31, 36] and in this sense Eq.(41) 
effectively takes the place of Schmidt decomposition. When 
|<7i) = |0) and \q 2 ) = |1), Eq. (42) gives the known answer for 
generalized GHZ state [25, 34]. 

The entanglement eigenvalue is minimal A^ ax = 1 /2 on 
condition that ki = k 2 . These states can be described as fol- 
lows 

|V> - |00 9l ) + \Uq 2 ) (43) 

where q\ and q 2 are arbitrary single qubit normalized states. 
The entanglement eigenvalue is constant A^ ax = 1/2 and 
does not depend on single qubit state parameters. Hence one 
may expect that all these states can be applied for teleporta- 
tion and superdense coding. It would be interesting to check 
whether this assumption is correct or not. 

It turns out that GHZ state is not a unique state and is one of 
two -parametric LU inequivalent states that have A^ ax = 1/2. 
On the other hand W-state is unique up to LU transforma- 
tions and the low bound Aj 2 nax = 4/9 is reached if and only 
if a = b = c. However, one cannot make such conclusions 
in general. Five real parameters are necessary to parameter- 
ize the set of inequivalent three qubit pure states [14]. And 
there is no explicit argument that W-state is not just one of LU 
inequivalent states that have A^ lax = 4/9. 

V. SUMMARY. 

We have derived algebraic equations defining geometric 
measure of three qubit pure states. These equations have a de- 



gree higher than four and explicit solutions for general cases 
cannot be derived analytically. However, the explicit expres- 
sions are not important. Remember that explicit expressions 
for the algebraic equations of degree three and four have a 
limited practical significance but the equations itself are more 
important. This is especially true for equations of higher de- 
gree; main results can be derived from the equations rather 
than from the expressions of their roots. 

Eq.(7) give the nearest separable state directly and this sep- 
arable states have useful applications. In order to construct 
an entanglement witness, for example, the crucial point lies 
in finding the nearest separable state [37]. This will be espe- 
cially interesting for highly entangled states that have a whole 
set of nearest separable states and allow to construct a set of 
entanglement witnesses. 

The expression in r.h.s. of Eq.(5) can be maximized directly 
for various three qubit states. Although it is very hard to solve 
the higher-degree equation, it turns out that the wide range of 
the three-qubit states have a symmetry and this symmetry re- 
duces the equations of degree six to the quadratic equations. 
In this reason Eq.(5) can be used to derive the analytic ex- 
pressions of the various entanglement measures for the three- 
qubit states. Also Eq.(5) can be a starting point to explore 
the numerical computation of the entanglement measures for 
the higher-qubit systems. We would like to discuss this issue 
elsewhere. 
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